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Fully discrete needlet approximation on the sphere^ 
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School of Mathematics and Statistics, UNSW Australia, Sydney, NSW, 2052, Australia 


Abstract 

Spherical needlets are highly localized radial polynomials on the sphere C d >2, with centers 

at the nodes of a suitable cubature rule. The original semidiscrete spherical needlet approximation of 
Narcowich, Petrushev and Ward is not computable, in that the needlet coefficients depend on inner 
product integrals. In this work we approximate these integrals by a second quadrature rule with an 
appropriate degree of precision, to construct a fully discrete needlet approximation. We prove that 
the resulting approximation is equivalent to filtered hyperinterpolation, that is to a filtered Fourier- 
Laplace series partial sum with inner products replaced by appropriate cubature sums. It follows 
that the Lp-error of discrete needlet approximation of order J for 1 < p < oo and s > djp has for 
a function / in the Sobolev space Wp(S'^) the optimal rate of convergence in the sense of optimal 
recovery, namely O (2“*^®). Moreover, this is achieved with a filter function that is of smoothness 
class c\-~ -I, in contrast to the usually assumed C°°. A numerical experiment for a class of functions 
in known Sobolev smoothness classes gives L 2 errors for the fully discrete needlet approximation that 
are almost identical to those for the original semidiscrete needlet approximation. Another experiment 
uses needlets over the whole sphere for the lower levels together with high-level needlets with centers 
restricted to a local region. The resulting errors are reduced in the local region away from the boundary, 
indicating that local refinement in special regions is a promising strategy. 
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1. Introduction 

Spherical wavelets [13] have wide applications in areas such as signal processing [20], geography 
[12, 35, 36] and cosmology [21, 33, 40]. The classical continuous wavelets represent a complicated 
function by projecting it onto different levels of a decomposition of the L 2 function space on the 
sphere. A projection, often called “a detail of the function”, becomes small rapidly as the level 
increases. This multilevel decomposition proves very useful in solving many problems. 

Narcowich et al. in recent work [29, 30] showed that the details of the spherical wavelets may be 
further broken up into still finer details, which are highly localized in space. This new decomposition 
of a spherical function is said to be a needlet decomposition. 
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Needlet approximation in its original form is however not suitable for direct implementation as its 
needlet coefficients are integrals. In this paper, we introduce a discrete spherical needlet approximation 
scheme by using spherical quadrature rules to approximate the inner product integrals and establish 
its approximation error for functions in Sobolev spaces on the sphere. Numerical experiments are 
carried out for this fully discrete version of the spherical needlet approximation. 

Before we describe spherical needlets and the discrete spherical needlet approximation we need 
some definitions. For d > 2, let be the real {d + l)-dimensional Euclidean space with inner 

product X • y for x, y G and Euclidean norm |x| := ^x • x. Let := {x G : |x| = 1} 

denote the unit sphere of The sphere forms a compact metric space, with the metric being 

the geodesic distance dist(x, y) := arccos(x • y) for x, y G 

Eor 1 < p < oo let Lp(§'^) = Lp(S^,(T(i), endowed with the norm |MlLp(§<i), be the real Lp-function 
space on with respect to the normalised Riemann surface measure ad on S'^. Eor p = 2, L 2 (S'^) 
forms a Hilbert space with inner product (/,5 )l 2(S‘^) /(^)5(^) dcTd(x) for f,gG L 2 (S'^). 

A spherical harmonic of degree £ on is the restriction to of a homogeneous and harmonic 
polynomial of total degree £ defined on Let T-Lf denote the set of all spherical harmonics of exact 

degree £ on The dimension of the linear space Idf is 


Z{d, £) 


{2£ + d-l) 


T{£ + d-l) 

r(d)r(£ + i) 


{£ + lf-\ 


( 1 ) 


where r(-) denotes the gamma function and ai x hi means cbi < ai < d hi for some positive constants 
c, c', and the asymptotic estimate uses [9, Eq. 5.11.12]. The linear span of ^ = 0,1,..., forms 
the space of spherical polynomials of degree up to u. Let be the Jacobi polynomial of 

degree £ for a, [3 > —1. We denote the normalised Legendre or Gegenbauer polynomial by 


p[d+l) 

-'T 


it) := Pi 


— ) ( — 



( 2 ) 


Given N > 1, for A: = 1,..., let x^ be nodes on and let Wk > 0 he corresponding 
weights. The set {{wk,Xk) : k = 1,... ,N} is a positive quadrature (numerical integration) rule exact 
for polynomials of degree up to u for some z/ > 0 if 


N 

/ p(x) dcrd(x) = y^u;fcp(xfc), 

ti 


for all p G 


Spherical needlets [29, 30] are a type of localized polynomial on the sphere associated with a 
quadrature rule and a filter. Let M+ := [0, +oo). 


Definition 1.1. A continuous compactly supported funetion g : M+ —)■ M+ is said to be a filter. We 
will only consider a filter with support a subset of [0, 2]. 

A filtered kernel on with filter g is, for T G M+, 


VT,gi^ • y) 


4,g(x-y) := < 


1 , 


e=o 


(fi+l) / 


'^ 9 ( 7 ^) Z{d,£)P) + Rx-y), 


0 < T < 1, 

T > 1. 


( 3 ) 
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We may define a filtered approximation on Li(S'^), T > 0 as an integral operator with the 
filtered kernel • y): for / G Li(S'^), 

^r,g(/;x) := l^T,g(/;x) := (/, r;r,g(x • •))L 2 (sd) = /(y) • y) d(Td(y). (4) 

Note that for T < 1 this is just the integral of /. 

We are now ready to define a needlet, following [29] who used a C°°(M+) filter and [30]. Let the 
needlet filter h be a filter with specified smoothness k > 1 (see Figure 1 in Section 5 for an example 
with K = 5) satisfying 


/i G C"^(M+), supp h = [1/2, 2]; (5a) 

h{tf+ h{2tf = 1 if te [1/2,1]. (5b) 


Condition (5b) is equivalent, given (5a), to the following partition of unity property for li^, 



-) =1. 

21) 


t > 1. 


For j = 0,1,..., we define the (spherical) needlet quadrature 


{{wjk,:sijk) ■ k = 1,.. Wjk>0, k = l,...,Nj, 

exact for polynomials of degree up to 2i~^^ — 1. 

A (spherical) needlet -ipjk, k = 1,Nj of order j with needlet filter h and needlet quadrature (6) is 
then defined by 

V'jfc(x) := U 2 i-i ;j(x • Xjfc), (7a) 

or equivalently, 

oo . 

'0ofc(x) := vVofe, '0jfc(x) := y) ^(d,^)L’y'^^^(x-Xjfc) if j > 1. (7b) 

t=o 

From (5a) we see that fijk is a polynomial of degree 2-^ — 1. It is a band-limited polynomial, in that 
'ipjk is L 2 -orthogonal to all polynomials of degree < 2-^“^. 

For / G L 2 (S'^), the original (spherical) needlet approximation with filter h and needlet quadrature 
(6) is defined (see [29]) by 

J N, 

x) := X] X] V'ifc(x), X G (8) 

j=0 k=l 

Note that Vj®®^(/;x) is a polynomial of degree at most 2*^ — 1 since ipjk is a polynomial of degree 
2i — 1. We shall call (/> V’ifc)L 2 (S‘*) semidiscrete (spherical) needlet coefficient and Vj‘®®‘^(/; •) the 
semidiscrete (spherical) needlet approximation to distinguish them from their fully discrete equivalents 
which we shall now introduce. 

The discrete (spherical) needlet approximation is defined by discretizing the inner-product integral 

/(y) V’ifc(y) do-d(y) 
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with another quadrature rule. For v > 0 and > 1, the discretization quadrature rule is 


Q^:=Q{u,N):={{Wi,yi):i = l,...,N}, Wi > 0, i = I,..., N, 
exact for polynomials of degree up to u. 

Let C'(S'^) be the space of continuous functions on S^. For f,g G C'(S'^), given Qat we define the 
discrete inner product by 

N 

{f,9)QN -=^^1 f{yi) g{yi). 
i=l 

The discrete (spherical) needlet coefficient of / for Qjv and i(jk is {f,'<Pjk)QN- We then define the 
discrete needlet approximation of order J by 

J 

:= xeS"*. (10) 

j=0 k=l 

Let h be a needlet filter with k > and let Wp(S^) with s > 0 and 1 < p < oo be a Sobolev 

space embedded in Lp(S'^), see Section 2.3 for the definition. In Theorem 4.3, we prove as a special 
case that for i/ = 3 • 2*^“^ — 1, i.e. Qn = Q{N, 3 • 2'^“^ — 1), the hp error using the approximation (10) 
for / G Wp(S‘^) and s > d/p has the convergence order i.e. 

11/- rj;r(/)|lL,(s-, < c2-'" ll/llw;(s-). / ew;(S'<). 

where the constant c depends only on d, s, h and n. This is consistent with the corresponding result 
for semidiscrete needlet approximation, see [29] and Theorem 3.12: 


||y _ <c2 /1| WKS'') • 

In Sections 3 and 4 we establish the connection to wavelets, and prove that the needlet approxima¬ 
tion is equivalent to a filtered approximation and that the discrete needlet approximation is equivalent 
to filtered hyperinterpolation [38] — a fully discrete version of the filtered approximation. These con¬ 
nections draw attention to the fact that the discrete needlet approximation considered in the present 
paper is not of itself new: what we have done is to express the filtered hyperinterpolation approxi¬ 
mation in terms of a frame {V’jfc} of the polynomial space where the frame has strong localization 
properties. The benefit will become apparent, however, if we take advantage of the local nature of 
the approximation to carry out local refinement. We make a preliminary study of local refinement of 
this kind in a numerical experiment in Section 5, though in this paper we do not develop the local 
theory. Rather, our main emphasis in this paper is on establishing the necessary theoretical tools 
for the discrete needlet approximation, on demonstrating the precise relationship between the various 
approximations, and on obtaining a global error analysis for / in Sobolev spaces. 

We note that Mhaskar [26, 27] proposed a full-discrete filtered polynomial approximation which 
is equivalent to filtered hyperinterpolation. A central assumption in [26, 27], in addition to polyno¬ 
mial exactness, is that a Marcinkiewicz-Zygmund (M-Z) inequality is satisfied. Quadrature rules with 
positive weights and polynomial exactness automatically satisfy an M-Z inequality (see Dai [7, Theo¬ 
rem 2.1] and Mhaskar [27, Theorem 3.3]). However, neither decomposition of wavelets into needlets 
nor numerical implementation were studied in [26, 27]. 
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Semidiscrete needlets are studied in other spaces, such as diffusion spaces and Dirichlet spaces, 
see [6, 17, 19, 24]. Quadrature rules and M-Z inequalities for these spaces are studied in [10, 11]. 
Establishing fully discrete needlet approximations on these spaces could be of great interest and will 
have wide applications. 

The paper is organised as follows. Section 2 gives necessary preparations. Section 3 studies the 
semidiscrete needlet approximation and its Lp approximation errors for / in Sobolev spaces on S'^, 
and its connection with the filtered approximation and continuous wavelets. In Section 4, we discuss 
the fully discrete needlet approximation and prove its approximation error for / G Wp(§'^) and exploit 
its relation to the filtered hyperinterpolation approximation and discrete wavelets. In Section 5, we 
give numerical examples of needlets and then some numerical experiments. Section 6 gives the proofs 
for the results in Section 3. 

2. Preliminaries 

The big O notation a{L) = Oa {b{L)) means that there exists a constant Cq, > 0 depending 
only on a such that \a{L)\ < Ca\b{L)\. The ceiling function [x] is the smallest integer at least x 
and the floor function [xj is the largest integer at most x. For integer k > 0 and real a > k, let 
(^) := = r(a-fc)r(^fc+i) extended binomial coefficient. We use “T” as a non¬ 

negative real number and use “L” as a power of 2, specially L := 2'^“^ for some J > 1, to remove any 
ambiguity about the meaning of the symbol L. 

2 .1. Lp space on sphere 

Let Lp(S‘^) = Lp(S'^,(Td) be the Lp-function space with respect to the normalised Riemann surface 
measure a = ad on endowed with the Lp norm 

ll/llLR§d) := |^J/(x)rdad(x)| , /GLp(S'^), 1 < p < oo; 

II/IIl^(s<^) := sup l/(x)i, / G Loo(s") nc(S''). 

xeS'* 

Since each pair of 77^, Hf, for 7 > 7' > 0 is L 2 -orthogonal, it follows that P,y(S'^) is the direct sum 
of T-if, i.e. = 0^=0Moreover, the direct sum 0^o'^r ^s dense in Lp(S'^), see e.g. [42, 

Ch.l]. Each member of Hf is an eigenfunction of the negative Laplace-Beltrami operator —A* on the 
sphere with eigenvalue 

= Af^=7(7 + d-1). (11) 

A zonal function is a function K : x —)• M that depends only on the inner product of the 
arguments, i.e. iL(x, y) = g((x • y), x, y G S'^, for some function g : [—1,1] —>■ M. Let • y) be 

the zonal function given by (2). From [39, Theorem 7.32.1, p. 168], 

|Ef+'^(x-y)| <1. (12) 

Let {Yi^m = '^i m : ^ > Oj m = 1,..., Z{d,i)} denote an orthonormal basis of the space L 2 (S'^). 
The normalised Legendre polynomial ■ y) satisfies the addition theorem 

Z{d,i) 

ye,m{^)Ye,m{y) = ^(d,QPf+'^(x • y). (13) 

m=l 
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This with the orthogonality of together gives the orthogonality of 

Let u(x • y) and ^(x • y) be two zonal functions of the form 


and 

t = P, 


oo oo 

v{y^-y) := Z{d,i) P^'^^^\yi-y), g{x-y) := Z{d,£) ■ y). 

£=o e=o 

Then (14) gives, for x, z G E>^, 

p OO 

(u(x • •),5'(z • •))l 2 (S‘^) = / • y) ^(z • y) daaiy) = J2^^be Z{d,£) z). 

JS‘‘ 


(14) 


(15) 


2.2. Sobolev spaces on the sphere 
Let s G M+. We define 

6^^^= (1 + A,)*/2 ^ (1 + ^)^ (16) 

where is given by (11). For 1" > 0, m = 1,..., Z(d, t), let 


flm ■— (/j ^,>ti)l 2 (§‘*) ■— / d(T£;(x) 


Sd 


be the Fourier coefficients of / G Li(S‘^). 

The generalised Sobolev space Wp(S'^) with s > 0 may be defined as the set of all functions 
/ G Lp(S'^) satisfying Ylm=P £ Lp(S'^). The Sobolev space Wp(S'^) is a Banach 

space with norm 


ll/llwj(§'^) 


OO Z{d,i) 

£=0 m=l 


Lp(§d)‘ 


(17) 


We have the following two embedding lemmas for Wp(S'^), see [18] and also [14, Eq. 14, p. 420]. Given 
a nonnegative integer k, let G''(S'^) denote the set of all k times continuously differentiable functions 
on 8“^. 


Lemma 2.1 (Continuous embedding into C"^(S'^)). Let d > 2. Given a non-negative integer k, the 
Sobolev space Wp(S'^) is continuously embedded into C'‘(S'^) if s > k + d/p. 

Lemma 2.2. Let d >2. For 0 < s < s' < oo and I < p < p' < oo, W®/(S'^) is continuously embedded 
into Wp(S‘^). 


2.3. Reproducing kernel Hilbert spaces 

When p = 2 and s > d/2 the space Wp(S'^) becomes a reproducing kernel Hilbert space. For brevity, 
we write BI''(§'^) := W|(S'^). The inner product in EI*(S'^) is (/, g)Hs(§d) := Ylm=P Tim dim- 

By (17) and L 2 -orthogonality of the II/IIhtS'^) = (E^o Em=f ^ • 

For s > d/2, each EI^(S'^) has associated with it a unique kernel K^^\x,y) satisfying, for x G 
and / G ]HI^(S'^), see e.g. [3, Section 2.4], 

iL('*)(x,.) GH^(S''), (/,7LW(x,.))^^^g^^ = /(x). (18) 
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From (18), we have 




x,y G 


(19) 


The kernel iF(^^(x, y) is said to be the reproducing kernel for EI^(S'^). It is moreover a zonal kernel, 
taking the explicit form 


(x • y) := V Z{d, i) pf(x • y) = 


£=0 


oo Z{d,£) 

EE- 

(=0 m=l 




Yi^rni^) Ye,miy)- 


( 20 ) 


3. Filtered operators, needlets and wavelets 

In this section, we study the properties of the filtered kernel, needlets and wavelets, and their 
relationships. 


3.1. Semidiscrete needlets and continuous wavelets 

We now point out the relation between spherical needlets and spherical wavelet decompositions. 
Let h he a needlet filter satisfying (5). Obviously the semidiscrete needlet approximation (8) can be 
written, for / G Li(§'^) and x G S'^, as 

,7 

yr'i(/;x) = j;W,(/;x), 

j=0 

where llj{f) is the contribution to the semidiscrete needlet approximation for level j: 

^7(/;x) := if,'^jkY(sd) V'ifc(x). (21) 

k=l 

In the language of wavelets, we may consider Uj(f]x.) to be the level-j “detail” of the approximation 

Needlets have a close relation to filtered polynomial approximations. At the heart of this relation¬ 
ship is the following expression, due to [29], and stated formally in Theorem 3.9 below: if ijjjk denotes 
the needlets of order j > 0 with needlet filter h and needlet quadrature (6), then 

Y'^jki^) i’jkiy) = ^ 2,-1 ;, 2 (x • y), (22) 

k=l 

in which the filter on the right-hand side, it should be noted, is /i^, the square of the needlet filter. 
This means that the level-j contribution to the semidiscrete needlet approximation can be written, 
using (21), as 

^i(/;x):= / /(x)u 2 y_i;, 2 (x-y) d(Trf(y). 

Js<^ 

To obtain the full semidiscrete needlet approximation we need to sum over j. For this purpose we 
introduce a new filter H related to the needlet filter h: 


H{t) : = 


1, 0<t<l, 

1 t > 1, 


(23) 
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and use the properties 


suppZ/C [0,2], =^h(^^y, t>l, JGZ+, (24) 

j=0 

which are easy consequences of (5). We note that this implies H G given h G C"'(M+). It 

then follows that X]/=o ^ 2^-1 ' y) = '^ 2 -^-^ ' y)^ J = 0,1,, and as a result the semidiscrete 

needlet approximation can be expressed as 

x) = / /(x) U 2 J -1 ^(x • y) dcJd(y) 

Js'i 

and for J > 1 this equals h(/)> where L := 2'^“^. 

3.2. Filtered operator and its kernel 

Recall the definition (3) of a filtered kernel. The convolution of two filtered kernels is also a filtered 
kernel. In particular, we have 

Proposition 3.1. Let d>2 and let g be a filter. Then for T > 0 and x, z G 

(^T,s(x • ■),VT,gi2 ■ •))L 2 (Sd) '= ' Y) ^T^gi^ ' Y) do-d(y) = -^^,92 ' z)- (25) 

Proof. For 0 < T < I, by (3), both sides of (25) equal I. We now prove (25) for T > 1. By (3) and 

(15), 

( CXD ^ CXD pf \ 

Z{d,i) Pf+^)(x ■■),Y,9{^) Z{d,i') • •) 

£=o e=o /L 2 (S'*) 

oo p 2 

= Z{d,e) P^‘^^'^\x-z) = Vj.g 2 {x-z), 

e=o 

thus completing the proof. □ 

When the filter is sufficiently smooth, the filtered kernel is strongly localized. This is shown in 
the following theorem proved by Narcowich et al. [30, Theorem 3.5, p. 584]. For integer k > 0, let 
C"^(M+) be the set of all k times continuously differentiable functions on M+. 

Theorem 3.2 ([30]). Let d > 2 and let g be a filter in C'‘(]R+) with 1 < k < oo sueh that g{t) is 
eonstant on [0, a] for some 0 < o < 2. Then 

|^t, 9 (cos 0 )| < T>1, (26) 

where cos 9 = x • y for some x, y G and the constant c depends only on d, g and k. 

We give an alternative proof of Theorem 3.2 in Section 6 , using different techniques. 

Remark. Dai and Xu [8, Lemma 2.6.7, p. 48] proved (26) for g G C'^'‘'’'^(M+). Brown and Dai [4, 
Eq. 3.5, p. 409] and Narcowich et al. ]29, Theorem 2.2, p. 533] proved that for g G C°°(M+), (26) 
holds for all positive integers k. 
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The following theorem shows the boundedness of the Li-norm of the hltered kernel. 

Theorem 3.3. Let d > 2 and let g be a filter in C"^(M+) with k > |_^^J such that g{t) is constant 
on [0, a] for some 0 < a < 2. Then 

lhTg(x • •)||li(s<^) - X e T > 0. (27) 

We give the proof of Theorem 3.3 in Section 6. 

Remark. Narcowich et al. [30] proved (27) for k> d + 1. 

Applying the convolution inequality of [2, Eq. 1.14, p. 207-208] to (4) gives 

II^T9(/)|lLp(Sd) ^ II^T,g(x- •)||Ll(§<^)ll/llLJ,(S‘^)• 

ThuS by Theorem 3.3, the operator norm of the filtered approximation on Lp(S'^) is bounded for 
g satisfying the condition of Theorem 3.2: 

Corollary 3.4. Let d > 2 and let g be a filter in C"'(M+) with k > |_^^J such that g{t) is constant 
on [0, a] for some 0 < a < 2 and let 1 < p < oo. Then the filtered approximation Vipg on Lp(S‘^) is an 
operator of strong type {p,p), i.e. 




'^^9 11 ILri— 


< Cd. 


.g,K, 


r > 0. 


(28) 


For L G Z+, the Lp error of best approximation of order L for / G Lp(S'^) is dehned by EL{f)p := 
^L{f)hp{sfi ■= 11/ “7'IIlp(S'*)' 

For given / G Li(S'^) and p G [l,oo], Fi(/)p is a non-increasing sequence. Since 
dense in Lp(S'^), the error of best approximation converges to zero as L —)> oo, i.e. limi_>.oo EL{f)p = 0, 
for / G Lp(S'^). 

The error of best approximation for functions in a Sobolev space has the following upper bound, 
see [18] and also [28, p. 1662]. 


Lemma 3.5 ([18, 28]). Let d > 2, s > 0 and 1 < p < oo. For L > I and f G Wp(S'^), 


EL{f)p<cL * ll/llw*(S'*)> 
where the constant c depends only on d, p and s. 


The filtered approximation Vp pj has a near-best approximation error for sufficiently smooth H in 
the sense of being within a constant factor of a best approximation error, as shown by the following 
lemma. 


Theorem 3.6. Let d > 2 and 1 < p < oo and let H be the filter given by (23) with h G C"^(M+) and 
K > |_^^J. Then for f G Lp(S'^) and L > 1, 

\\f-VL,Hif)LpiS<^)^^^L{f)p, (29) 

where the constant c depends only on d, H and k. 
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Proof. Since H{t) = 1 for t G [0,1], we have 

•))L2(sd) =f^(x), 

and hence, from Corollary 3.4, 

||/-^L,//(/)||lj,(S'^) = \\f ~P~^L,Hif-p)\\hpisd) < (1+||^L,//|Ilp^Lp)II-^“^'IIlp(S‘*) ^ Cd,H\\f-p\\hpiS<i)^ 

which holds for all p G Pi(S^), thus completing the proof. □ 

Remark. The estimate (29) is a generalisation of the results of Rustamov [34, Lemma 3.1, p. 316] 
and Sloan [37]. Rustamov proved (29) for H G C'°°(M+) and 1 < p < oo while Sloan showed (29) for 
p = oo, f G C(S'^) and H G C"^'''^(M+), and even for certain piecewise polynomial filters H belonging 
to C‘^-i(M+). 

Lemma 3.5 and Theorem 3.6 give the error of the filtered approximation for Sobolev spaces: 
Corollary 3.7. With the assumptions of Theorem 3.6, for f G Wp(S'^) with s > 0 and L > 1, 

11/ - <cL ^ ll/llw-(Sd), 

where the constant c depends only on d, p, s, H and n. 

Remark. H. Wang and Sloan in a recent work [fl] proved Corollary 3.7 using essentially the same 
method. 


3.3. Semidiscrete needlet approximation 

The smoothness of the filter makes the needlet ^[jk localized. This can be seen from the following 
corollary of Theorem 3.2, first proved by Narcowich et al. in [30, Corollary 5.3, p. 592]. 


Corollary 3.8 ([30]). Let h be a needlet filter, satisfying (5). If h ^ C'^(M+) with k> d + 1, then 

c2J'^ 


IV’jfcWI < 


, X G S'^, j > 0, A: = 1,.. .,Nj, 


(1 + 2l dist(x,Xjfc))” 
where the constant c depends only on d, h and k. 


The following theorem shows, as foreshadowed in (22), that an appropriate sum of products of 
needlets is exactly a filtered kernel. It is implicit in [29]. 

Theorem 3.9 (Needlets and filtered kernel). Let h be a needlet filter, see (5), and let H be given by 
(23). For j > 0 and 1 <k < Nj, let i[ijk be needlets with filter h and needlet quadrature (6). Then, 

^ V’ifc(x) V'ifc(y) = n 2 i-i ;j 2 (x • y), j>0, (30a) 

k=l 
J Nj 

EE V’jfc(x) i’jkiy) = • y), J > 0. (30b) 

j=0 k=l 

For completeness we give a proof. 
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Proof. For j = 0, by (7a) and (3) 


No No „ 

V V'ofc(x) ^ofc(y) = y' = / dcjrf(z) = 1 = t> 2 -i ^2 (x • y). 

For j > 1, using (7b) and the fact that the filter h has support [1/2, 2], we have (noting h{2) = 0) 

Nj 2J-12J-1 „ 

^V'ifc(x) V’,fc(y) = 

fc=l £=0 £'=0 

X ^Wjk Z{d,£)P^/^^\x ■ Xjk) Z{d,£')P^f"^^\y ■ xjk). (31) 

k=l 

Since {(rCjfc,Xjfc) : k = l,...,Nj} is exact for polynomials of degree 2-^+^ — 1, the sum Ylk=i 
quadrature points in (31) is equal to the integral over Then by (14) and the definition of the 

filtered kernel, see (3), the equation (31) gives 


oo p 2 

X^V’ifc(x) V’ifc(y) = ^id,£) = ^2,-1 ;j2(x • y). 

k=l £=0 


This proves (30a). 

For J > 0, by (32) and (24), we now have, using /i(0) = 0, 


J 

EE ipjk{x) V’ifc(y) 

j=0 k=l 


J 00 p 2 

l + Z{d,£)p'f^^\i^-y) 

j=i i=i 

‘+e('E'<I) 1 z(<!,f)pr''(x-y) 

i=l \j=0 J 

00 p 

e=o 


This completes the proof. 


(32) 


□ 


Theorem 3.9 with (4) leads to the following equivalence of the filtered approximation with filter 
H and the semidiscrete needlet approximation (8). 

Theorem 3.10. Under the assumption of Theorem 3.9, for f G Li(S'^) and J > 0, 

J 

= E E (/• '/»)l,(S*) ’/’Jl = ><?'"“(/). (33) 

i=o k=i 

Theorems 3.6 and 3.10 imply that the semidiscrete needlet approximation has a near-best approx¬ 
imation error. 

Theorem 3.11. Let J G Z+ and let V)/®®‘^(/), see (8), be the semidiscrete needlet approximation with 
needlets V’jfc, see (7), for filter smoothness n > |_y^J. Then for 1 < p < 00 and f G Lp(S'^), 

where the constant c depends only on d, the filter h and k. 
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Proof. By Theorem 3.10, the approximation by the semidiscrete needlets is equivalent to that 

by filtered approximation Then the definition (8) of and (29) of Theorem 3.6 

together with Theorem 3.10 give 

11/ - K?'"‘(/)||l,(S-) = 11/ - ''2^-..«(/)|Il,(S-) S E22-.(/),, 

□ 

Theorem 3.11 and Lemma 3.5 imply a rate of convergence of the approximation error of Vj®®'^(/) 
for / in a Sobolev space, as follows. 

Theorem 3.12. Under the assumptions of Theorem 3.11, we have for f G Wp(S‘^) with s > 0 and 

J > 0, 

I|/-K,”"'(/)|Il,(s-|<<^2-'-||/||„.„.,, 

where the eonstant c depends only on d, p, s, h and k. 

4. Discrete needlet approximation 

To implement the needlet approximation in a numerical computation, we need to discretize the 
continuous inner product (/, in (8). We make use of the quadrature rule in (9) to replace the 

continuous inner product by a discrete version. In this section, we estimate the error by the discrete 
needlet approximation for the Sobolev space Wp(S‘^), 1 < p < oo. 

4 . 1 . Discrete needlets and filtered hyperinterpolation 

Let be needlets satisfying (7), and let Qn := Q{N,i) := {{Wi,yi) : i = 1,...,A^} be a 
discretization quadrature rule that is exact for polynomials of degree up to some I, yet to be fixed. 
Applying the quadrature rule Qn to the needlet coefficient (/, V’ifc)L 2 (§d) = /s'* /(y)V*j 7 c(y) d(Td(y), we 
obtain the discrete needlet coefficient 

N 

if,'f>jk)QN fiyi) i^jkiyi)- (34) 

i=l 

This turns the semidiscrete needlet approximation (8) into the (fully) discrete needlet approximation: 

J 

^jk)QN (35) 

j=0 k=l 

In a similar way to the semidiscrete case, cf. (33) of Theorem 3.10, the discrete needlet approxi¬ 
mation (35) is equivalent to filtered hyperinterpolation, which we now introduce. 

The filtered hyperinterpolation approximation with a filtered kernel Vj< g in (3) and discretization 
quadrature Qat in (9) is 

N 

VT,g,Nif',^) := ^T,g,iv(/;x) := {f,VTg{- :='^Wif{yi)vTg{yi-y), T G M+, (36) 

i=l 

as named by Sloan and Womersley [38]; see also [22] and [16]. 
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Theorem 4.1. Let h he a needlet filter given by (5) and let the filter H he given by (23) . For f € C (S*^) 
and J > 0, 

J N, 

^2^-1,/^,7v(/) = = VjXif)- (37) 

j=0 k=l 

Remark. Note that in Theorem fil we do not yet require the number N of nodes of the discretization 
quadrature to depend on the order J of the discrete needlet approximation. 

Proof. Applying (30b) of Theorem 3.9 to ^ 2^-1 //(y* ' (36), and using (34), we have 

N J Nj 

i=l j=0k=l 

J Nj / N \ J Nj 

= EE E Wi f(yi) ipjkiyi) j f’jk(x) = EE(/ ) '4^jk')Qiq fi’jki^)^ 
j=0 k=l Vi=l / j=0 k=l 

which gives (37). □ 

4 . 2 . Error for filtered hyperinterpolation 

By Theorem 4.1, the discrete needlet approximation, if regarded as a function over the entire sphere, 
reduces to the filtered hyperinterpolation approximation. In this section, we deduce the approximation 
error of the filtered hyperinterpolation or discrete needlet approximation for / in Sobolev spaces 
Wp(S'^) with 1 < p < 00 and s > d/p. 

For comparison, by Corollary 3.7, the continuous approximation nif) following approx¬ 

imation error for / G Wp(S'^) with 1 < p < 00 and s > 0: 

11/“ ^l,_h'(/)IIlp(S‘*) — c 7/ ||/|lw®(S‘*)) /£Wp(S'^), (38) 

where the constant c depends only on d, s, filter H and k. 

We now want error bounds for For that discrete version of the filtered approximation, 

Mhaskar [27, Theorems 3.1 and 3.2] obtained the same truncation error to that of (38) for / G Wp(S'^) 
with 1 < p < 00 , see also [22, 38] which proved the case when p = 00 , and s > d/p, but with a stronger 
assumption on the filter smoothness than stated in Theorem 4.2 below. Given J G Nq, let L := 
and let 

Qn := Q{N, 3L - 1) := {(W,, y,) : i = 1,2,..., N} (39) 

be a discretization quadrature exact for polynomials of degree up to 3L — 1. The following theorem is 
due to H. Wang and Sloan [41]. 

Theorem 4.2. Let d>2, l<p<oo and s > d/p. Given a needlet filter h, let Vi^ ^ N 
filtered hyperinterpolation in (36) with Qn given by (39) and filter H given by (23) and satisfying 
H G G^(M+) for K > . Then, for f G W;(S'^), 

11/ - ^L,/^,Ar(/)|lLj,(S‘') < cT ^ ll/llw^(§‘i); 
where the constant c depends only on d, p, s, H and n. 
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Theorem 4.2 with Theorem 4.1 gives the errors for the discrete needlet approximation of / G 
Wp(S'^), 1 < p < oo, as follows. 

Theorem 4.3 (Error by discrete needlets for Wp(S'^)). Let d >2, 1 < p < oo and s > d/p, and let 
be the discrete needlet approximation given by (10) with needlet filter h G C"‘(M+) and k > |_^^J 
and with diseretization quadrature Qn in (39). Then, for f G Wp(S'^) and J G No, 

11/ - K7^?^'^(/)|Ilp(S'*) < c 2 ll/llwj(§‘^)) 

where the eonstant c depends only on d, p, s, h and k. 


4-3. Discrete needlets and discrete wavelets 

Let Qat be a discretization quadrature rule given by (9). The discrete needlet approximation 
in (10) can be written, for / G C'(S'^) and x G S'^, as 

J 

rj;U‘‘{/;x) = (40) 

j=0 

where is the level-j contribution of the discrete needlet approximation defined by 

Nj 

L^i, 7 v(/;x) :=L/j- 7 v.(/;x) := f G C(S'^), x G S'^. (41) 

k=l 

Using (30a) then gives 

= (/, J]]'0jfc(-)V’ifc(x))^^ = (^/,U2i-i ;,2(x- = ^2,-1 (42) 

k=l 

where the hltered kernel U 2 j_i ^ 2 (x • y) is given by (3). 

Using (30b) and (42) with (4) gives the following representation of filtered hyperinterpolation in 
terms of L(j,N- 


Theorem 4.4. Let d > 2 and let Uj^isf{f) be the level-j eontrihution of the discrete needlet approxi¬ 
mation in (40) and let H be the filter given by (23). Then for f G C'(S‘^) and J G No, 

J 

j=0 

Theorems 4.1 and 4.4 with (36) and (42) imply the following representation for 

Corollary 4.5. Let h be a needlet filter given by (5) and let the filter H be given by (23). For 
f G C(S'^) and J G No, 

J 

= ^ J]]lUi/(yi)u2,-q;,2(yi-x) = (/,U2j-q^(-•x))g^. (43) 

j=0 i=l 

The theorem below shows that the Lp-norm of decays to zero exponentially with respect 

to order j. This means that the different levels of a discrete needlet approximation have different 
contributions and Uj^M{f) thus forms a multilevel decomposition. We can hence regard Uj^Nif) as a 
discrete wavelet transform. 
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Theorem 4.6. Let d>2, l<p<oo and s > d/p and let be the level-j contribution of the 
discrete needlet approximation in (40) and let the needlet filter h satisfy h G C"^(M_|_) and k > |_^^J • 
Then for f G Wp(§'^) and j > I, 

where the constant c depends only on d, p, s, h and k. 

Proof. Theorem 4.4 shows that Uj^N{f) is the difference of two filtered hyperinterpolation approxima¬ 
tions: for j > 1, 

This with Corollary 4.2 gives 

< c2 ||/||w|(S'*)) 

where the constant c depends only on d, p, s, h and n. □ 

5. Numerical examples 

In this section we give a computational strategy for discrete needlet approximation and show the 
results of some numerical experiments. For the semidiscrete needlet case the approximation is not 
computable, but we are able to infer the error indirectly by using the Fourier-Laplace series of the 
test function to evaluate the L 2 error. The last part gives an example of a localized discrete needlet 
approximation with high accuracy over a local region. 

5.1. Algorithm 

Algorithm 5.1. Consider computing the discrete needlet approximation Vj^‘^(/;x') of order J G No 
with needlet filter h at a set of points {x' : i = 1,... ,M}. The needlet quadrature rules {{wjk,:>i.jk) ■ 
k = l,...,Nj} are exact for polynomials of degree 2^^^ —1. The major steps are analysis and synthesis. 

1 . Analysis: Compute the discrete needlet coefficients (/, 'fijk)QN , fe = 1,..., Nj, j = 0,..., J using a 
discretization quadrature rule Qn = Q(A, 3 • 2*^“^ — 1) = {(VFj, y^) : i = 1,..., N}. 

2. Synthesis: Compute the discrete needlet approximation i = 1, • • •, M. 

Needlet filters. Here fijki'idf) is computed by (7b) where the normalised Legendre polynomial 
p^(d+i)(^) computed by the three-term recurrence formula, see [9, § 18.9(i)] and the needlet fil¬ 
ter may be computed as follows. For construction of other needlet filters, see e.g. [25, 29]. 

Given k > 1, let p{t) be a polynomial of degree 2 k -|- 2 of the form 

2k+2 

p{t)-= X] ^e[0,1], (44) 

k=K-\-l 

where the coefficients are uniquely determined real numbers satisfying p(0) = 1 and the ith deriva¬ 
tives of p{t) at t = 0 for 1 < f < K -|- 1 are zero. Clearly, p[l) = 0 and all the jth derivatives of p(t). 
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Figure 1: Needlet filter h e (7^(11+) Figure 2: An order-6 needlet with a C^-needlet filter 


ait = l are zero. Then it can be shown that 

h{t) := < 


l<t<2, 


P{t - 1 ), 

0, elsewhere 


is a hlter h satisfying (5). This section uses k = 5, where the coefficients in (44) are: ag = 924, 
ay = —4752, 08 = 10395, og = —12320, oio = 8316, on = —3024, 012 = 462, giving the filter h 
illustrated in Figure 1. 

Figure 2 shows an order-6 needlet with the filter given in Figure 1. We see that it is very localized. 
Quadrature rules. We use symmetric spherical designs for integration on S^, as recently developed 
by Womersley [44], for both the needlet quadrature rule and the discretization quadrature rule. Let t 
be a non-negative integer. A symmetric (if Xj is a node so is —Xj) spherical t-design is a quadrature 
rule with equal weights and exact for all polynomials of degree at most t. In these experiments the 


t-designs have 2 


t'^+t+4 


t^/2 nodes so are more efficient than t-designs using [t + 1)^ (which is 
ynomial space Pt(S^)) points, see for example [1]. The symmetric spherical 


the dimension of the po 

designs also have good geometric properties [44], which will beneht localized needlet approximations 
(see Section 5.3 below), compared with other quadrature rules such as [15] on §^. 

Cost of algorithm. Using a symmetric spherical t-design, a needlet quadrature rule for level j has 
Nj nodes, giving a total of X]/=o ~ I ^ 2^'^ nodes for all J levels as the symmetric spherical 

t-designs are not nested. Similarly, a discretization quadrature rule exact up to degree 3 x 2'^~^ — 1 
has Af ~ I X 2^'^ nodes. Thus the analysis step to evaluate the needlet coefficients requires | x 2^"^N 
evaluations of /. The synthesis step only involves a weighted sum of the needlets evaluated at M 
(possibly very large) points. At high levels the number of needlets is large, for example when J = 6, 
L = 64, Nj = 8130 and N = 4562. 


5.2. Needlet approximation for the entire sphere 

This section illustrates the discrete needlet approximation of a function / that is a linear combi¬ 
nation of scaled Wendland radial basis functions on S^, see [43]. The advantage of this choice is that 
the Wendland functions have varying smoothness, and belong to known Sobolev spaces. 
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Let (r)+ := max{r, 0} for r G M. The original Wendland functions are [43] 


( 1 - 0 ^, 

k = 0, 

(l-r)^(4r + l), 

k = 1, 

(l-r)®(35r2 + 18r + 3)/3. 

k = 2, 

(1 - r)%{32r^ + 25r2 + 8 r + 1 ), 

k = 3, 


(1 - r)]^°(429r^ + 450r^ + 2Wr‘^ + 50r + 5)/5, k = 4. 


The normalised (equal area) Wendland functions as defined in [5] are 


■ — 4^k 



:= 


(3fc + 3)r(fc + i) 
2r{k + l) 


k>0. 


The Wendland functions scaled this way have the property of converging pointwise to a Gaussian as 
k —)> oo, see Chernih et al. [5]. Thus as k increases the main change is to the smoothness of /. We 
write 4 >(r) := for brevity if no confusion arises. 

Let zi := (1,0,0), Z 2 := (-1,0,0), Z 3 := (0,1,0), Z 4 := (0,-1,0), Z 5 := (0,0,1), zq := (0,0,-1) 
be six points on and define [23] 


6 

/(x) :=/fc(x) := ^(/>fc(|zi - x|), A: > 0, 
i=l 


(45) 


where | • | is the Euclidean distance. 

Narcowich and Ward [31] and Le Gia, Sloan and Wendland [23] proved that fk G (S^). 
Figure 3 shows the picture of / 2 , which belongs to B[ 2 (S^). The function fk has limited smoothness 
at the centers Zj and at the boundary of each cap with center Zj. These features make fk relatively 
difficult to approximate in these regions, especially for small k. 

L 2 approximation error. We show the L 2 errors when using and by For we 

compute its L 2 error by discretizing the squared L 2 -norm by a quadrature rule. We cannot compute 
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the L 2 error for in this way as we do not have access to exact integrals for the inner products. 

As the test function in (45) is a linear combination of Wendland functions, we are able to approximate 
the L 2 error of by truncating the Fourier-Laplace expansion and using the known Fourier 

coefficients of Wendland functions. 

We make use of the Fourier-Laplace coefficients of / to compute the L 2 -error of the semidiscrete 
needlet approximation over the entire sphere, as follows. By Theorem 3.10 and the definition of the 
filtered approximation, see (3) and (4), and the addition theorem, see (13), the Fourier coefficients of 
yjieed(j) fgjn- Then the Parseval’s identity gives 

00 2£-|-l 

II/-k;‘~‘(/)||Ls=,= e E {1 -»(46) 

f=2-4-l+l m=l 


We expand 0(^2 — 2t) in terms of Pi{t): 

00 

4>{V2^t) = Y,^g {2e+l)Pe{t), 
e=o 

where Pi{t) is the Legendre polynomial of degree i and 

cj,g:=-j ^ 4>{V2^t)Pg{t) dt, £ > 0. (47) 

Using the addition theorem again, 

00 00 

(/)(|zi -x\) = (j) (V2^^Yz7^) = (2^ + P)Pi{'2‘i ■ x) = EE 

£=0 e=0 m=l 


which with (45) gives 


6 

fim if 

i=l 


This with (46) and the addition theorem together gives 


1/ - Kf'"'{/)|lLs-) 


00 2£dr\ 

E E(i-«(At)) 

£=2-4-1+1 m=l 


£ \\^ \2 |2 


6 6 


^ ^ F£,m (2 


i=l 


|<)),|"j;j;(2£ + l)P,(z*-z,), 

e=2J-l+l i=l j=l 


(48) 


where we use the Gauss-Legendre rule to compute the one-dimensional integral (47) for 4£ to the 
desired accuracy. 

Figure 4a shows the L 2 -error of the semidiscrete needlet approximation for A: = 0,1, 2 ,3,4, 

where we used the filter h of Figure 1, with H then given by (24), and the order of semidiscrete needlet 
approximation is J = 1,..., 6, and the truncation degree i in (48) is taken as high as 500. The slight 
fluctuation of the L 2 -errors of the semidiscrete needlet approximation for is partly due to the trun¬ 
cation error for the Fourier coefficients of 4a- We use a log-log plot with the horizontal axis indexed 
by orders J and fit data from J = 3 to 6 for each k to illustrate the convergence order. 
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10° 



V k = 0, s = 1.5 
- 7.0e-01 2 '’'-^'’ 

* k = 1,s = 2.5 
- 9.7e+00 2'®'® 

* k = 2, s = 3.5 
2.1e+02 2'®®® 

« k = 3, s = 4.5 

- 7.7e+03 2'^®® 

O k = 4, s = 5.5 
- 4.06+05 2'^““® 



Order J 

(a) Semidiscrete (b) Fully discrete 

Figure 4: L 2 -errors for needlet approximations of fk, with orders J = 1, ..., 6, filter h G 


Either (43) and (3), or the needlet decomposition (10) can be used to compute the fully discrete 
needlet approximation Some discussion of efficient implementation can be found in [16]. We 

then approximate the L 2 error by a quadrature rule {{wi, Xj) : i = 1 ,..., as follows. 

N 

- /(x)f dcjrf(x) - /(Xi))^ (49) 

Figure 4b shows the corresponding L 2 -error for the discrete needlet approximation 
where we used the same needlet filter, and used symmetric spherical designs for both needlets and 
discretization, and the order of discrete needlet approximation is J = 1,... ,6. We used a symmetric 
spherical 301-design (with N = 45454 nodes and equal weights Wi = 1/N) to approximate the integral 
in (49). 

For each k, the L 2 -errors of the semidiscrete and fully discrete needlet approximations converge 
at almost the same order (with respect to 2'^). For the particular test function fk the order of 
convergence is close to better than the optimal order when considering all functions in 

EI^(S^). The figure also shows that the convergence order becomes higher as the smoothness of / 
increases, which is consistent with the theory. 

5.3. Local approximation by discrete needlets 

In the following example, we show the approximation error using discrete needlets for /2 given by 
(45), using all needlets at low levels and needlets with centers in a small region at high levels. 

In general, let X be a compact set of S'^. We define the localized discrete needlet approximation 
for / E C(S^) by 

Jo 

xeS'^Vx, 

j=0 
Jo 

(/>V’ifc)sjvV’jfc(x), xeX, 

j=0 
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Figure 5: Centers of needlets at level 4 (larger 
points) and level 6 (smaller points) for a local¬ 
ized discrete needlet approximation 



1 10 


8 


6 


4 


2 


I 


Figure 6: Absolute errors of a localized dis¬ 
crete needlet approximation for /2 with full 
needlets at levels < 4 and localized needlets at 
levels 5,6 


where i^j_Ar(/; x), given by (41), is the level-j contribution of the discrete needlet approximation. 
The idea is that on the compact set X we seek a more rehned needlet approximation — that is, we 
“zoom-in” on the set X. 

Let X := C (z 3 ,r), the spherical cap with center Z 3 := (0,1,0) and radius r. 

Figure 6 shows the pointwise absolute error of the localized discrete needlet approximation 
V"/‘g®^(C (z 3 , 7 r/ 6 ) ;/ 2 ;x). In the exterior of the cap C (z 3 , 7 r/ 6 ), the approximation used needlets up 
to level 4, with the largest absolute error, about 1.2 x 10“^, at the centers Zj, i 7 ^ 3. In the cap, the 
approximation is a combination of needlets at the low levels 0 to 4 with those at high levels 5 and 6 . 

We observe that the localized discrete needlet approximation has good approximation near the 
center of the local region but with less computational cost since the levels 5 and 6 used only a fraction 
of the full set of needlets, approximately \C (Z 3 , r) |/|S^| = (1 — cos(r ))/2 (about 6.7% when r = vr/G). 
This localization is an efficient way of constructing a discrete needlet approximation for a specific 
region. 

Figure 5 shows the centers of the needlets for level 4 (larger points) and those in the cap C (Z 3 , tt/G) 
for level 6 (smaller points) of the localized discrete needlet approximation ’l 4 "g®^(C (Z 3 , vr/G);/ 2 ; x). 
The smaller points illustrate where the high levels of the localized discrete needlet approximation 
focused. At level 6 , the needlet quadrature used the symmetric spherical 63-design, which has totally 
8130 nodes over the sphere and 544 nodes in the cap. The localized discrete needlet approximation at 
this level used only needlets with centers at these 544 nodes for the local region. 
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6. Proofs for Section 3 


In the proof of Theorem 3.2 we use the following lemma to bound Mi^(T,£), where 


Mk{T,i) : = 


9 


9 


Mk-iiT,i) Mk-i(T,e + l) 


k = 1, 

/c = 2,3,.... 


( 50 ) 


2 i + 2 r + k 2 {i+l) + 2 r + k' 

Note that vanishes for f < \aT~\ — k, because of the assumed constancy of g on [0,a]. This 

is a crucial property for establishing the following lemma. 

Lemma 6.1. Let g satisfy the condition of Theorem 3.2 with Ti = \aT~\ and T 2 = [2TJ and with T 
sufficiently large that 0 < Ti — n < T 2 . Let Mk{T, i) be defined by (50). Then for an arbitrary positive 
integer k < k, 

Mk{T, e) = 0 , Ti-k<i<T 2 , (51) 

where the constant in the big O depends only on d, k, g and k. 

Proof. For a sequence ui, let Aj U£ := := U£ — U£+i denote the first order forward differ¬ 

ence of U£, and for i > 2 , let the ith order forward difference be defined recursively by A\{u£) := 
a] (A^“^(u£)). We now prove the estimate in (51), making use of the obvious identity 

'A] {U£ V£) = (^^ U£) V£ + U£+l ()^] V£). (52) 

By (50), for A; > 2 


Mk{T,l) = 


Mk-iiT,£) _ 

2 i + 2r + k 2{£ + l) + 2r + k 


Mk-iiT,£) 


+ 


Mk-i{T,e) 


Mfc i(r,f + i) 


2 {i + l) + 2r + k 2{i + l) + 2r + k 


+ 1} ) Mk-i{T,i) =: 6k{i){Mk-i{T,i)). 


2 i + 2r + k + 2 \2£ + 2r + k 
In addition, let Si{i) := A]. Then for /c > 1, 

By induction using (52) and (53), with k > 1 can be written as 


Mk{T,£) = Y,R-i2k-i-i){e) Alg(^-)., 

i=l 


(53) 


(54) 


where A: — 1 < j < 2A: — 2, is a rational function of A with degree^ deg(i?_j) < — j and hence 


R-j{t) = Od,k [r^] 


(55) 


For g G C'^(M+) and 0 < i < k < k, we have by induction the following integral representation of 
the finite difference A\g{^)-. 


AW 


T 


dui 


. j_ 
T 


,(*) 


-I- Ui H- \-Ui dui. 


^Let R{t) be a rational polynomial taking the form R{t) = p{t)/q{t), where p{t) and q{t) are polynomials with q fi 0. 
The degree of R{t) is deg(i?) ;= deg(p) — deg(g). 
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Since is bounded, for Ti — k < i < T 2 , 




< Ci,g 


This together with (54) and (55) gives (51), on noting that £ x T in (51). 


(56) 

□ 


Proof of Theorem 3.2. In this proof, let r := (d — 2)/2, Ti := \aT] and T 2 := [2rj. We only need to 
consider T sufficiently large to ensure that 0 < Ti — k < r 2 . Let t G [—1,1], be the Jacobi 

polynomial of degree i for a,f3 > —1. From [39, Eq. 4.5.3, p. 71], 

[2j + a + p + i) L\j + a + p + i) ^ 


E 

j =0 


(2j + a +l))r(j + a + 13 + l) p(a,/3)^^^ _ r(£ + a + /3 + 2) 


r(j + /3 + i) r(£ + /3 + i) 

and by [39, Eq. 4.1.1, p. 58], pj“’^^(l) = (^”^“). Then we find using (1) and (2) that 

OQ p 

VT^gicosO) = Pi‘^^^\cose) 


(57) 


e=o 


ni) 

rid) 


d\ 00 


E^ 

e=o 


^ \ L^{+ + 2.^ +il p{r,r) 


r(|) 

m 


Y^M^iTp 


Tie + r + 1) 

Ti^ + 2 r + K + 1 ) ^{r+K,r) 


£=Tl —K 


r(£ + r + l) 


(cos 0 ). 


(58) 


where the last equality uses (57) and summation by parts n times, and MfiiT,i) is given by (50). 
Erom [39, Eq. 7.32.5, Eq. 4.1.3] or [ 8 , Eq. B.1.7, p. 416], for arbitrary a, (3 > —1, 


p1“’^^(cos6 *)| < —— 
^ ' ('^-1 


^a,/3 ^ ^ 


0 < 9 < TT. 


1 ) ^ _ 


(£-i+ 0 )«+|(£-i+vr- 0 )^+i 
Applying Lemma 6.1 with (58) and (59) gives (bearing in mind that r = [d — 2)/2) 

\vTg[eose)\<Cd,^ V lM,(T,£)ir+" X - 

T 2 

E Cd q K / j_. j_. . 

(7-i + 0)-+xr(7-i + ^-0)^ 

From this and Ti x T x T 2 together with T 2 — Ti x T, for 6 G [0, 7 r/ 2 ] we have 

'J'd 

1U7^„(COS0)1 < CrfoK p-{2K.-l) - CcdqK -XT; 

' (T-i + 9)-+"-^ - (1 + T9)^+"-^ 


(59) 


(r 1 + 9 )^+'^+2 (7-1 + TT - 0)^+2 

p-{2K-l) 


while for 9 G [7r/2,7r], 


I '^T,g \ 


e=i i - 


The estimates for the above two cases imply (26), thus completing the proof. 
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Proof of Theorem 3.3. In the proof, let Ti := \aT] and T 2 := [2TJ. We only need to consider T 
sufficiently large to ensure that 0 < Ti — k < T 2 . Using summation by parts k times, for t G [—1,1], 


CXJ ^ ^‘2 / p \ 

£=0 £=Ti-k ^ ^ 

where for k G Z+, := where we used an asymptotic estimate of Gamma function, 

see e.g. [9, Eq. 5.11.13, Eq. 5.11.15], and 


1 


Pj i e [-1,1] 

j=o 

is the Kogbetliantz (Cesaro) kernel of order k, see e.g. [42, Section 2.3] and [2, Eq. 2.1.5]. Berens et 
al. [2, Theorem 2.1.1], see also [42, Theorem 2.3.10], proved that for k > and x G S'^, 


O' 


hv—l 


(x-OI 


Li(S‘^) — '^d.,K,g- 


This together with (60) and (56) gives 


Ta 




T 


1^/^—1 I II 1 




ILi(§^) 


i=Tl —K 

< x; r-r-n,,,., 

e=Ti-K 




thus completing the proof. □ 
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